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Abstract. Because of the spatially long-ranged nature of spontaneous fluctuations in thermal non¬ 
equilibrium systems, they are affected by boundary conditions for the fluctuating hydrodynamic variables. 
In this paper we consider a liquid mixture between two rigid and impervious plates with a stationary con¬ 
centration gradient resulting from a temperature gradient through the Soret effect. For liquid mixtures with 
large Lewis and Schmidt numbers, we are able to obtain explicit analytical expressions for the intensity of 
the non-equilibrium concentration fluctuations as a function of the frequency co and the wave number q of 
the fluctuations. In addition we elucidate the spatial dependence of the intensity of the non-equilibrium 
fluctuations responsible for a non-equilibrium Casimir effect. 

PACS. 05.40-a Fluctuation phenomena, random processes, noise, and Brownian motion - 05.70.Ln Ther¬ 
modynamics - Nonequilibrium and irreversible thermodynamics - 65.40.De Thermal properties of crys¬ 
talline solids - Thermal expansion; thermomechanical effects - 68.60.Dv Physical properties of thin films, 


nonelectronic - Thermal stability; thermal effects 

1 Introduction 

It has been well established, both theoretically and exper¬ 
imentally, that thermal fluctuations in fluids in the pres¬ 
ence of a temperature gradient and/or a concentration 
gradient are spectacularly long range wm- Specifically, 
the intensity of the temperature and concentration fluc¬ 
tuations in non-equilibrium (NE) fluids varies with the 
wave number q of the fluctuations as q ~ 4 which means 
that in real space the correlations extend over the size 
L of the system [J. As a consequence, the spatial spec¬ 
trum of these NE fluctuations is strongly affected by finite- 
size effects due to the presence of boundaries, especially 
for small wave numbers accessible in shadowgraph exper¬ 
iments. Hence, for the interpretation of NE shadowgraph 
experiments [mmciEiiaiiniiiiiiii] and for the validation 
of computational studies |13U14U151il61ll7j of NE fluctu¬ 
ations, an assessment of the impact of the boundary con¬ 
ditions of the fluctuating hydrodynamic variables on the 
NE fluctuations is necessary. 

A procedure for solving the fluctuating hydrodynam¬ 
ics equations to obtain the intensity of the NE concentra¬ 
tions fluctuations has been developed by two of us, but 
with artificial boundary conditions for the fluctuations at 
the walls adopted for mathematical convenience )18j . For 
the case of realistic boundary conditions, we subsequently 
obtained an approximate solution in terms of a Galerkin 
approximation Ca¬ 


in liquid mixtures there are two diffusion modes that 
are linear combinations of heat and mass diffusion and a 
viscous mode wm- Important parameters are the Lewis 
number, that is the ratio of the thermal diffusivity Dt and 
the mass diffusion coefficient D 1 and the Schmidt num¬ 
ber, that is the ration of the kinematic viscosity v and D. 
In liquid mixtures both Lewis and Schmidt numbers are 
commonly larger than unity, implying that temperature 
and viscous fluctuations decay much faster than concen¬ 
trations fluctuations. Hence, in dealing with liquids one 
often adopts an approximation of large Lewis and Schmidt 
numbers, Le 1 and Sc 1. This approach is particu¬ 
larly convenient when the focus is on concentration fluctu¬ 
ations at time scales when temperature and viscous fluc¬ 
tuations are fully decayed j51[2i?j. For Le 1 and Sc 1 
we were able to obtain an exact expression for the au¬ 
tocorrelation function of the intensity of the NE concen¬ 
tration fluctuations in the presence of realistic boundary 
conditions, but that required a tedious numerical evalua¬ 
tion of the decay rates (eigenvalues) of the modes of the 
hydrodynamic operator, yielding some analytical results 
only in the long-wavelength limit m- Subsequently, the 
procedure was extended to study the dynamics of the NE 
concentration fluctuations and good agreement with ex¬ 
perimental measurements was obtained, although again a 
tedious numerical evaluation of the appropriate eigenval¬ 
ues was required YU 1 . 
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NE concentration fluctuations are not only affected by 
finite-size effects, but also by the presence of gravity. In 
this paper, as in some earlier work [15], we focus our at¬ 
tention specifically on the finite-size effects. We have been 
able to obtain relatively compact analytic expressions for 
both the dynamic and static autocorrelation function of 
the NE concentration fluctuations that do require only 
minimal numerical work. Our results are not only rele¬ 
vant to past and future experimental studies of NE con¬ 
centration fluctuations at low gravity [T0lf25] , but also for 
certain aspects of earth-bound non-equilibrium phenom¬ 
ena, like a NE Casimir effect [2511371125] . as well as for the 
interpretation of numerical simulations of non-equilibrium 
fluctuating hydrodynamics in mixtures Emm]. 

We shall proceed as follows. The equations for the ap¬ 
propriate fluctuating variables and boundary conditions 
are specified in Section [2] The concentration gradient in¬ 
duces a coupling between concentration fluctuations and 
wall-normal velocity fluctuations along the direction of the 
concentration gradient mm- The equation for these 
velocity fluctuations is solved in Section [3] Inserting this 
solution in the equation for the concentration fluctuations 
we obtain the correlation functions for the NE concentra¬ 
tion fluctuations in Section [5] In Section [5] we analyze the 
spatial dependence of the intensity of the NE concentra¬ 
tion fluctuations needed for an understanding of Casimir 
forces induced by the NE concentration fluctuations [25 . 
Our results are summarized and discussed in Section [5] 


2 Fluctuating variables and boundary 
conditions 

We consider a layer of a binary fluid mixture bounded by 
two plane-parallel walls separated by a distance L and lo¬ 
cated at z = ± 7 jL. The plates are maintained at different 
temperatures, so that a stationary temperature gradient 
VT = (VT) z exists, with z the unit vector in the wall- 
normal direction. As a consequence of the Soret effect, 
a stationary concentration gradient Vc = (Vc) z will 
also appear, where c is the concentration in mass frac¬ 
tion of one of the components of the mixture. We shall 
further consider that all relevant thermophysical proper¬ 
ties of the mixture are constant (do not depend on tem¬ 
perature or concentration) so that both VT and Vc are 
uniform. This physical situations is usually referred to as 
the Rayleigh-Benard problem for a binary mixture 1221132] . 
In the absence of gravity, the stationary state described 
above is stable. Finally, consistent with the Le 1 ap¬ 
proximation, we assume the separation ratio of the mix¬ 
ture if) = —/ 3Vc/aVT to be positive [221155] (the most 
common case). Here a = — 1/p (dp/dT) c ^ p is the thermal 
expansion coefficient and (3=1 /p ( dp/dc)T,p the solu- 
tal expansion coefficient, with p the mass density of the 
mixture. 

Fluctuating hydrodynamics, originally developed for 
dealing with spontaneous thermodynamic fluctuations in 
equilibrium [5511551155] . can also account for fluctuations 
around the stationary gradient specified above [51155], In 


the most general case, as in equilibrium [55], there will 
be coupled pressure, temperature, concentration and ve¬ 
locity fluctuations, so that some simplifications are desir¬ 
able. If Sv are the velocity fluctuations, one first assumes, 
as usual for dense fluids, that they are divergence-free: 
V • Sv (incompressible flow). This assumption allows us 
to neglect the sound modes [30]. Next, we adopt here a 
large Lewis number approximation (Le 1) that decou¬ 
ples temperature and concentration fluctuations. Hence, 
one has to consider only the coupling of concentration and 
wall-normal velocity fluctuations 55112511551 . The Le > 1 
approximation was first proposed to simplify the linear 
stability analysis of the Rayleigh-Benard problem in bi¬ 
nary fluids mm, demonstrating its validity for positive 
separation ratios. Subsequently, it was used for simplifying 
the associated fluctuating-hydrodynamics problem mm 
133] , It has been successfully employed in numerical sim¬ 
ulations EE for reproducing light-scattering experiments 
in binary mixtures [24] and, more recently, in ternary mix¬ 
tures EM?]. We note that for the most common liquid 
mixtures, Le is of the order 10 — 10 3 . 

The fluctuating hydrodynamics equations for our prob¬ 
lem have been presented in some previous publications, to 
which we refer for details and physical background EH 
mm- For our purpose here we find it convenient to use 
the working equations of Ref. 25] in dimensionless form. 
Hence, we introduce dimensionless space, time and con¬ 
centration variables by: 


r = r/L, 



Sc = 


Sc 

Vc L' 


( 1 ) 


In terms of these dimensionless variables, and in the ab¬ 
sence of gravity, the fluctuating hydrodynamics equations 
for a binary fluid mixture in the Le 1 limit become [251 : 


-^d t (V 2 6v z ) = V 4 5v z -F 1 (r,t) 

DC 

d t Sc = V 2 Sc - Sv, z + T 2 (r, t) 
where Sc = v/D is the Schmidt number and 

Fi(r,t) = (V x V x (V • 5TI)} z , 


pvD 


F 2 {r 1 t) = - 


pDVc 


V -SJ. 


(2a) 

(2b) 


(3) 


We note that in Eqs. ©-©, for simplicity, tildes in di¬ 
mensionless variables have been suppressed. From here on, 
everything will be in terms of dimensionless variables, ex¬ 
cept where explicitly the opposite is stated. 

As in our previous work [23] , to further simplify Eqs. © 
we also assume the Schmidt number to be large, Sc 1. 
Hence, in practice, the left-hand side (LHS) of Eq. (I2al) 
will be taken as zero. As Le 1, Sc 1 is also a good 
approximation for typical liquid mixtures. This approxi¬ 
mation amounts to assuming that viscous fluctuations de¬ 
cay very fast in the time scale set by the diffusion time, 
see Eq. ©. 
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Equations m define dimensionless random forces ex¬ 
pressed as spatial derivatives of the two fluctuating dis¬ 
sipative fluxes: a random stress tensor 5II(r,t ) and a 
random diffusion flux SJ(r,t). We note that, in fluctuat¬ 
ing hydrodynamics [3J1331M1I3MSS1 , the linear phenomeno¬ 
logical relations, used to ’close’ the balance laws and to 
obtain the hydrodynamic equations, need to be supple¬ 
mented with fluctuating dissipative fluxes. In our case, 
a sn is added to Newton’s viscosity law, and a 6J is 
added to Fick’s diffusion law. The statistical properties of 
these fluctuating dissipative fluxes are given, in terms of 
the dissipation matrix, by the so-called Fluctuation Dis¬ 
sipation Theorem (FDT). Since the random diffusion flux 
contributes only to equilibrium fluctuations, in what fol¬ 
lows only the FDT for the random stress tensor will be 
needed, which for the case of incompressible flow and in 
terms of dimensional variables, it explicitly reads mm-. 

(^■(r,t) Sn* kl (r',t')) =2k B T V 

x {SikSji + SuSjk ) 5{t - t') 5(r - r'). (4) 

where k b is Boltzmann’s constant and 77 the shear vis¬ 
cosity. As previous investigators mm we assume that 
a local version of the FDT continues to be valid in non 
equilibrium, but we neglect the effects of nonhomogeneous 
noise so that, in practice, temperature T in Eq. 0 is iden¬ 
tified with the average temperature in the layer. 

Realistic boundary conditions are no-slip for the ve¬ 
locity fluctuations and impervious wall for the mass flow. 
For an incompressible fluid, and in terms of dimensionless 
variables, they are [401 : 

0 = 6v z = d z Sv z , 0 = d z Sc, at, z = (5) 

In the second of the boundary conditions above the limit 
Le ~^> 1 is implicit, meaning that temperature fluctua¬ 
tions are neglected and consequently the thermodiffusion 
contribution to the mass flow. 

Next, to solve Eqs. @ with the boundary conditions 0 , 
as usual one applies Fourier transforms in 

time and in the AW-plane, parallel to the walls. Hence, in 
matrix form and for Sc 1, Eqs. © become: 


\dl - q 2 ) 2 0 

Sv z 


Fi(u},q,z) 

1 i uj + q 2 — d 2 

6c 


F 2 (u,q 1 z) 


where to is the frequency of the fluctuations and q is the 
component of the wave vector of the fluctuations in the 
horizontal plane (with corresponding wave number q). The 
random forces Fi[u>, q, z) in Eq. © are the Fourier trans¬ 
forms of the random forces in Eq. ©. To solve Eq. © we 
first separate its right-hand side (RHS) as [57] : 


Fi{uJ,q,z) 


0 

I 

Fi(u;,g,z) 

F 2 {bJ,q,z) 


F 2 (u,q,z) 

I 

0 


This splits the solution of Eq. © additively into two parts: 
Sc = <5c E + 5c NE that we distinguish with superscripts ’E’ 
and ’NE’. It is interesting to note that the random stress 
tensor 5II and the random diffusion flux 5J are uncorre¬ 
lated, so that the random forces F± and F 2 in the RHS 


of Eq. © are uncorrelated too, see Eq. ©. As a conse¬ 
quence, the autocorrelation function of the concentration 
fluctuations splits additively also into ’E’ and ’NE’ com¬ 
ponents. Furthermore the ’E’ component of the autocor¬ 
relation is the same as for the system in equilibrium, i.e., 
for Vc = 0 (or the 1 were not present in the first column 
of the hydrodynamic matrix of Eq. ©). From now on we 
concentrate on the ’NE’ contribution to the concentration 
fluctuations to be determined from Eq. ©, namely: 


\(d 2 ~q 2 ) 2 0 

5v z 


Fi{u:,q,z) 

1 iuj T q 2 — d 2 

5c 


0 


where, for simplicity, we dropped superscript ’NE’. Of 
course, Eq. © must be solved subjected to the bound¬ 
ary conditions © in the wall-normal z-variable. 


3 Wall-normal velocity fluctuations 


We start by solving the velocity equation, the first of 
Eqs. ©. As in previous works, the solution is conveniently 
expressed in a series of eigenfunctions I4l fl9ll23ll39] 

OO 

6v z (z) =J2 Cn w *i z )i (9) 

JV=1 

where the hydrodynamic modes Wn(z) are the solutions 
of the eigenvalue problem 

{dl - q 2 ) 2 w N (z) = r 2 w N (z), (10) 

subject to the boundary conditions 

W N (z) = d z Wn{z) = 0, at z = ±\. (11) 

In Eqs. m-m we have used the fact that the eigenprob- 
lem mi has an infinite numerable set of solutions Wn{z) 
with distinct eigenvalues T^ that are real and positive 
numbers. Although we have explicitly calculated the so¬ 
lution to Eqs. CDD-CD in terms of hyperbolic functions, 
including the numerical solution, for each q , of an alge¬ 
braic equation to obtain r 2 r , we shall not present it here. 
It is very technical and, as clarified below, we only need 
the trace of the differential operator in the LHS of Eq. m 
with boundary conditions ED. As happens in many prob¬ 
lems in physics of fluids mm, the trace of the hydrody¬ 
namic operator can be computed without detailed discus¬ 
sion of the hydrodynamic modes. 

Indeed, independently of the details of Wn{z), it fol¬ 
lows that: 


, 1/2 

/ dzW M (z)[(d 2 -q 2 ) 2 W N (z)\ = 
J- 1/2 


, 1/2 

/ dz[(d 2 -q 2 ) 2 W M (z)\ W N {z ), (12) 

1 - 1/2 


obtained upon integration by parts and use of the bound¬ 
ary conditions ED- As a consequence, the Wn(z) form an 
ortogonal set, namely 


, 1/2 

/ dz Wm(z) Wn(z) = D n Snm- 
J- 1/2 


(13) 
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Hence, from Eqs. m and m the coefficients Cn in the 
series expansion m can be readily evaluated, resulting in: 


C N 


F ltN (u,q) 

D N (q) r$r(qy 


(14) 


with 


f 1/2 

Fi,jv(w,q) = / dz W N (q,z) Fi(w, q,z). (15) 

J- 1/2 

This procedure solves the first of Eqs. © with the appro¬ 
priate boundary conditions. 

However, we are interested only in the two-point cor¬ 
relation function ( Sv z (z ) 5v* z (z')). By using Eq. ([9]) it can 
be expressed as a double series of hydrodynamic modes 
Wn(z). From Eq. (HUH) it follows that the coefficients of 
this double series are related to the correlation matrix 
(Fi ; jv(w, q) Fj" M (u >’, q')}. This correlation matrix, in turn, 
can be calculated from the definition of Fi(r,t), Eq. ©, 
and the FDT for the random stress 511 (r, t), Eq. fl3}. This 
last calculation has been already presented in detail pre¬ 
viously (3ll42| . Here we simply quote the final result: 

(Fl,at(w, q) F* tM (uj',q')} = F NM (q) 

x (27 t) 3 5(oj — u/) 5(q — q'), (16) 


of the velocity fluctuations were retained in the LHS of 
Eq. Ball if 5] , the the wall-normal velocity autocorrelation 
would depend explicitly on ui. 

As in previous publications mm, the sum of the se¬ 
ries (1201) can be performed analytically. Indeed, applying 
the operator (d 2 — q 2 ) 2 at the two members of Eq. (BUI) 
one has: 


(d 2 z ~ q 2 ) 2 V(q , z, z') = 2 Fq 2 S(z - z'). (21) 

Notice that, since 5(z — z') as a function of z satisfies the 
boundary conditions m, it can be expanded in a series 
of eigenfunctions Wn(z), namely 


S(z-z') = J2 


N =0 


F>N{.q) 


W N (q,z)W N (q, Z '), (22) 


independent of q , from which Eq. m follows. Hence, if 
one solves Eq. BT1) subject to the boundary conditions 

V(q,z,z') = d z V(q,z,z')=0, at z = ±\, (23) 

which follows from Eq. ED. we obtain an analytic expres¬ 
sion for V(q, z, z') presented in AnnendixlAllsee Eq. (1A41I) ). 


4 Non-equilibrium concentration fluctuations 


with 


F N M{q) = 2,Fq 2 

= 2Fq 2 rjj{q) D N {q) Snm- 


1/2 

dzdz’ [D 2 W N (q,z )] ^(?',2'), 
- 1/2 


(17) 


where, in the second line of Eq. ED. we used Eqs. ED 
and (fl3l) . In Eq. 1171) V 2 = (d~ — q 2 ) 2 , while the symbol 
F denotes a dimensionless prefactor 


P = ML. 

pvDL 


(18) 


Now combining all the information above we are able to 
obtain an explicit expression for the wall-normal velocity 
autocorrelation function. In view of Eq. ED one has: 


To solve the second of Eqs. © for the concentration fluc¬ 
tuations with the appropriate boundary conditions, we fol¬ 
low an equivalent approach to the one in the previous sec¬ 
tion for the wall-normal velocity fluctuations. That is, we 
expand 5c in a series of eigenfunctions of the differential 
operator d 2 — q 2 subject to the boundary conditions. In 
this case this approach is quite simple since the eigenfunc¬ 
tions are trigonometric functions. Hence, one has: 


5c(uj,q,z) 


-A 0 ^-~v A n cqs{2Nitz) 

i u) + q 2 iuj + q 2 + AN 2 tt 2 

_ 9 V- Bn sin[(27V + 1)7rz] 

^ iw + q 2 + (2N + 1) 2 tt 2 1 j 

AT=0 v ’ 


with 


(Sv z (u],q,z) 5v* z (w',q',z')) = V{q, Z> z') 

x (27 t) 3 5(co — w') 5(q — q r ), (19) 


with 

°° 2Fo 2 

V{q,z,z') = n ( \ s W N{q,z) W N (q,z’), (20) 
D N (q) r^{q) 

where we use the fact that the eigenvalues and eigenfunc¬ 
tions are real and positive. It is interesting to note that 
the autocorrelation V(q, z, z') in Eq. ED c fo es n °t depend 
explicitly on the frequency u. This is a consequence of the 
Sc 1 approximation, which neglects the inertial term 
in the starting equations, that is, the LHS of Eq. (Ball 
is assumed to be zero [231124] . If the temporal derivative 


A N {u;,q) 
Bn (w, q) 


,1/2 

I - 1/2 
,1/2 


dz 5v z (uj, q, z) cos(2Nttz), 


,i/Z 

/ dz 5v z (ui,q, z) sin[(2IV + l)7rz]. 

J-1/2 


(25) 


From Eq. HMD we see that the autocorrelation function 
(5c(uj 1 q 1 z)5c*(u} , ,q' 1 z')) of the NE concentration fluctu¬ 
ations is directly related to (5v z (uj,q,z) 5vt(oj', q\ z')), 
the autocorrelation function of the wall-normal velocity 
fluctuations discussed in the previous section. In view of 
Eq. ED ; the result can be cast in the form 


(6c(u>,q,z) 5c* (u/, q', z')) = S{u,q,z,z') 

x (27 t) 3 5(uj — uj') 5(q — q'), (26) 
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where S(ui, q , z, z') is a double series of trigonometric func¬ 
tions. For generic z and z the explicit expression for 
S(u), q, z, z') is quite lengthy, so that we present it explic¬ 
itly in Appendix [B] see in particular Eq. (1B44I) . Fortu¬ 
nately, for practical applications, one does not need the 
general S(uj, q, z, z'), but some particular values and/or 
integrals that allow for more compact expressions. 

In light-scattering [5115T115T] or shadowgraph mm ex¬ 
periments one commonly studies the NE fluctuations with 
the wave vector in the horizontal direction. The struc¬ 
ture factor S(ui,q) that is measured in such experiments 
is given by a double integration over the thickness of the 
layer of the S(u, q, z, z') of Eq. (OBI) , namely: 



Integrating Eq. (1B44I) for S(oj, q, z, z 1 ), one readily ob¬ 
serves that only the first term of the series (that does 
not depend on z or z') gives a non-zero contribution to 
S(ui, q). Hence, one simply has: 


S(v,q) 


2 F 1 4(1 — cosh q) 

to 2 + q 4 \_ q 2 q 3 (q + sinh q) 


(28) 


Applying double inverse Fourier transforms to Eq. (1261) . 
and taking into account Eq. 423), one obtains the inten¬ 
sity S/g) of the equal-time NE concentration fluctuations 
observed in experiments: 



= S{q)(2nf 5(q-q'). (29) 

Integration of Eq. (1551) over the frequency oj yields a com¬ 
pact exact analytic expression for 5(g) 


S(q) = 


2tt 


= F 


doj S(ui, q) 

! 

4(1 — cosh < 7 ) 


g 4 g 5 (g + sinhg) 


(30) 


From Eq. HU we obtain at large q 1 


S{q)^ 



(31) 


whose first term contains the typical 5(g) ~ g -4 behav¬ 
ior of NE temperature and concentration fluctuations dis¬ 
cussed in many publications [1 13(11144] . For this particular 
problem of NE concentration fluctuations induced by the 
Soret effect, the g -4 behavior was first found by Law and 
Nieuwoudt [3D] without gravity and by Segre and Sen- 
gers [U] with gravity. In the limit q —> 0, 5(g) from 
Eq. m reaches a finite limit, namely 


in agreement with Ref. m for a solutal Rayleigh number 
equal to zero (no gravity). 

All results presented so far are in terms of the dimen¬ 
sionless variables defined in Eq. ©■ In terms of dimen¬ 
sional variables our primary result, the NE structure fac¬ 
tor given by Eq. (1551) . becomes: 

1 4(1 — coshg) 

g(g + sinh g) 

. (33) 

where g = qL. In obtaining Eq. (1331) . Eq. 43 f° r dimen¬ 
sionless concentration and Eq (1151) for F, were used. 

It is interesting to observe, either in Eq. (1281) or in 
Eq. (1331) . that in the absence of gravity the NE concen¬ 
tration fluctuations have a single decay time, even when 
boundary conditions are included in the calculation. This 
is in contrast with the situation when gravity is not ne¬ 
glected, for which a series of exponentials is obtained for 
the NE time correlation function M- The physical rea¬ 
son behind this difference is that gravity affects the decay 
time of NE concentration fluctuations, not only the ampli¬ 
tude, as first discussed by Segre and Sengers m- That in 
microgravity NE concentration fluctuations have a single 
decay time, even at very small g, as predicted by Eqs. (1551) 
or (l33l) , has indeed been observed experimentally |9]. 

The exact result of Eq. hu supersedes a Galerkin ap¬ 
proximation proposed some years ago |19] which, in the 
limit of zero gravity, was used for the analysis of micro- 
gravity experiments mm- Comparison of the two ex¬ 
pressions shows that both share the same asymptotic limit 
for large g, while the Galerkin approximation [19] overes¬ 
timates the exact small g limit of Eq. (1551) by ~ 30%. 


5(w,g) = 


k B T (Vc ) 4 2Dg 2 
p uDq A ui 2 + D 2 q 4 


5 Correlations relevant for the NE Casimir 
effect 

Because of their spatially long-ranged nature, it has been 
recently proposed that NE fluctuations will induce forces 
(Casimir forces) somewhat similar to the Casimir forces 
appearing at equilibrium critical points. To evaluate these 
NE Casimir forces [2611271155] one expands the pressure up 
to second order in the fluctuating fields, so that the result¬ 
ing fluctuation-induced force becomes proportional to the 
mean square intensity of the NE fluctuations in real space. 
For the case of a binary mixture [551 , the mean square in¬ 
tensity of concentration fluctuations that is relevant to 
Casimit forces, can be related to the structure factors dis¬ 
cussed in this paper by first considering S(q,z), obtained 
from the two-points S(uj,q,z,z r ) of Eq. (1561) . upon inte¬ 
gration over the frequency to and substitution of z = z', 
namely 

1 

5 '( 9 , 2 ) = —J duj 5(w, q, z, z). (34) 

This quantity is related to the equal-time mean square NE 
concentration fluctuation at a single spatial point: 


S(q) 


F 

720 


Fq 2 

15120 


+ 0(q 4 ) 


(32) 


(6c(r, f) 2 ) N E = ( 1 5c(~) 2 ) 


^2 J 2 ?rg S(q, z) dq , (35) 
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Fig. 1. Intensity of the NE concentration fluctuations S(q, z), 
relative to the value of the structure factor S(q), as a function 
of 2 for three wave numbers q. In the limit q —> 0, S(q, z ) and 
S(q) become identical. The deviations increase with increasing 
q, but since the contributions for large q fall off as l/q 4 , S(q) 
is a reasonable first approximation in practice. 

obtained by applying inverse Fourier transforms in the two 
frequencies and in the two (2D) wave vectors to Eq. (1251) . 
Notice that (<5c(r, £) 2 )ne is both stationary (it does not de¬ 
pend on t) and translationally invariant in the XU-plane, 
i.e., in a direction parallel to the bounding walls (it does 
not depend on the coordinates x , y of the point r). It only 
explicitly depends on the coordinate 2 of the point r. 

From Eq. (IB44I) for q, z, z'), an explicit expression 
for S(q 7 z) can be readily obtained, namely: 

.A A nm cos(27Vttz) cos(2Mnz) 
S(g,*) = SM + 2££ 

N=0 M—l * 

Bnm sin[(2iV + 1)7 tz \ sin[(2M + 1)7 rz] 

+ ^ > 2 

Q 2 + \ [(2iV + l) 2 + (2 M + l) 2 ] 

(36) 

with matrices A^M^q) and B^M^q) defined by Eqs. (1B45I) 
and (IB46I) . respectively, and S(q) being the ‘experimental’ 
static structure factor given by Eq. m - Equation (1561) 
for S(q , z) is a long expression that can only be simplified 
marginally. Hence, we continue our discussion of S(q,z) 
graphically. 

Figure [T] shows the intensity of the NE concentration 
fluctuations S(q,z) of Eq. (1551) . relative to the value of 



Fig. 2. NE ( Sc(z ) 2 ) as a function of 2 evaluated from Eqs. (1351) - 
(l36l) (red line) together with NE (5T(z) 2 ) for rigid boundaries 
(blue line) as evaluated in previous work [27] for Pr = 6. For 
an easy comparison both curves are normalized independently, 
so that the average value in the layer is, for each case, equal to 
unity. Note that ( ST(z ) 2 ) vanishes at the walls, while ( 5c(z ) 2 } 
does not. The arrows represent the approximation (|37[l for 
(5c 2 ) wa n. 


the structure factor S(q) of Eq. (fBOl) . as a function of 
z for three wave numbers q , as indicated. One first ob¬ 
serves that, just as for temperature fluctuations in a one- 
component fluid [27], S(q,z) has an inhomogeneous spa¬ 
tial distribution that reaches a maximum at mid-layer. 
However, we stress that, in contrast to the case of NE 
temperature fluctuations, S(q,z ) for the NE concentra¬ 
tion fluctuations does not vanish at the walls, 2 = ±4- 
The reason is that to establish a uniform temperature 
gradient, in experiments the fluid layer is confined be¬ 
tween thermally conducting plates EHEUEQL so that the 
temperature fluctuations, in contrast to the concentration 
fluctuations, must vanish at the walls of the plates mm- 

Upon substitution of Eq. (1M1) into Eq. (1551) one can 
obtain an explicit expression for the intensity of NE con¬ 
centration fluctuations (c>c( 2 ) 2 ). The resulting integrals in 
q do not admit a compact analytical expression, so that 
we also discuss this quantity graphically. In Fig. [2] the red 
curve represents the quantity {5c{z) 2 ) computed by eval¬ 
uating numerically the q -integral in Eq. (1351) . It is inter¬ 
esting to compare it with the intensity ( 5T(z ) 2 ) of the NE 
temperature fluctuations in a one-component fluid that 
was investigated in an earlier publication m for realistic 
rigid boundary conditions. Actually, ( ST(z ) 2 ) depends on 
the Prandl number Pr [451 : the blue curve in Fig. [2] shows 
{5T(z ) 2 ) for Pr = 6 [27| . For easier comparison we normal¬ 
ized the two quantities independently, so that the average 
value through the layer (area below the curve) is, in each 
case, equal to unity. Figure [2] shows the main conclusion of 
this section: The different boundary condition (vanishing 
gradient vs. vanishing function) causes different behaviors 
at the walls, NE temperature fluctuations vanish while 
NE concentration fluctuations do not. We also emphasize 
that (6c(z) 2 ) reaches a non-zero value at the walls with a 
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non-zero slope, although maybe not completely evident in 
Fig. [2] because, due to numerical limitations, only a finite 
number of modes was added. 

To finalize, it is interesting to evaluate the layer-average 
square of concentration fluctuations, that is the quantity 
actually relevant for Casimir forces [2Ej. Upon substitu¬ 
tion of Eq. (Idiil) into Eq. (Idl^l one obtains: 


f r°° 

1 4(1 — cosh q) 

2tt Jq q 

q 4 q 5 (q + sinhg) 


= 3.11 x 1(T 3 F, 

in terms of dimensionless variables. To revert the result to 
physical dimensions one has to use Eq. m for dimension¬ 
less concentration and Eq m for F. so that: 

(Sc 2 ) = 3.11 x 10~ 3 -^L(Vc) 2 . (38) 

puD 

We have indicated with arrows the result m in Fig. [3 
with the same normalization as ( Sc(z ) 2 ). One observes a 
small difference between the actual value of (i5c 2 ) wa ii and 
the value given by Eqs . EB-®. 

6 Summary and conclusions 

In this paper we have investigated NE concentration fluc¬ 
tuations in a binary fluid mixture in the absence of gravity. 
We adopted a realistic setting where a stationary concen¬ 
tration gradient is induced, through the Soret effect, by 
maintaining the two bounding plates of a fluid layer at 
different temperatures. No-slip boundary conditions were 
used for the velocity fluctuations and no-flow for the con¬ 
centration fluctuations. An approximation of large Lewis 
and Schmidt numbers was adopted, which means that 
both temperature and velocity fluctuations decay much 
more rapidly than the concentration fluctuations, so that 
the slow dominant mode is a pure concentration mode. 
This is a good approximation for dense liquid mixtures, 
while more questionable for rarefied gases. 

With all the features described in the previous para¬ 
graph, we have been able to obtain exact relatively com¬ 
pact analytical expressions for both the dynamic structure 
factor, Eq. ®, and the static structure factor, Eq. ®, 
that would be observed in heterodyne low-angle light scat¬ 
tering or shadowgraph experiments. These expressions will 
be useful for the interpretation of experimental results, in 
particular in microgravity conditions as well as for com¬ 
puter simulations of NE concentration fluctuations m- 
We completed the paper studying the total intensity of 
the NE concentration fluctuations as a function of the dis¬ 
tance to the walls. We demonstrated the important result 
that, due to the different boundary condition, the inten¬ 
sity of concentration fluctuations does not vanish at the 
walls, in contrast to the case of temperature fluctuations. 
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A Explicit expression for the wall-normal velocity autocorrelation function 


A particular solution of Eq. m is: 


V p (q,z,z') 



z'\ cosh(g |2 


- sinh(g|z 



(A39) 


The general solution of the homogeneous problem corresponding to Eq. m (t-e., by taking the RHS equal to zero) 
is: 

Vh(q, z, z') = A 0 cosh qz + Ai sinh qz + A 2 z cosh qz + A 3 Z sinh qz, (A40) 

where the coefficients At will be, in general, arbitrary functions of z' and q. Adding Eqs. (1A39H and (1A40I) . and 
imposing the four boundary conditions (E^l) . the coefficients Ai(q, z') are uniquely determined. We thus obtain an 
explicit expression for V{q,z,z'), namely 


V(q,z,z') = 


where 


F 

L 2 

, 1 a 2 

q 2 — sinh 2 q 

l 

4 q 3 . 


q sinh q 


A i (q) = 


zz 


cosh[g(z — z')\ + A 2 (z — z') sinh[g( 2 : — z')\ 

A! 


4 q 2 sinh q 
sinh q + q cosh q 


cosh[g( 2 i + z')\ - Ai(z - z') sinh[g(z - /)] j + V p (q,z,z'), (A41) 


A 2 (q) = 


q + sinh q cosh q 


(A42) 


We have checked, by numerically evaluating the decay rates Ay and the modes Wn(z) solving Eq. (flOl) . for two 
different values of q, that indeed Eq. (1A41I) is equal to the sum of the series EU1) . 


B Explicit expression for the concentration autocorrelation function 

Here, we explicitly present the double trigonometric series for the two-point autocorrelation function S(ui,q, z, z') of 
the NE concentration fluctuations, introduced in Eq. (ini). We first note that there is a slight simplification, since for 
the function V(q, z, z') of Eq. (1A41I) : 

, 1/2 , 1/2 

/ dz dz' V(q,z,z') cos(2Nttz) sin[(2M + 1)ttz\ = 0 (B43) 

J- 1/2 J- 1/2 

for any pair of integers N, M. Thus, S(u, q 1 z, z') reads explicitly: 


S(cv,q,z,z') = 


2 F 


1 4(1 — cosh q) 


lo 2 + q 4 [q 2 q 3 (q + sinhg) 


2Aiyo{q) cos(2Nwz) 

^ (—iw + q 2 ){ito + q 2 + AN 2 / k 2 ) 




+ A) TV, 


2A]y 0 (q) cos(2Nttz') 


T.Y.k 


4Ajvm(?) cos(2 Nnz) cos(2Mnz') 


— ' (ico + q 2 )(—iuj + q 2 + 4:N 2 n 2 ) 2 -—' ' (iw + q 2 + AN 2/ n 2 )(— iw + q 2 + 4M 2 7r 2 ) 
jv=i v ^ ^ ' JV=1M=1 \ 1 1 / 

4sin[(27V + 1)7 tz\ sin[(2M + 1)7 tz'] 


OO OO 


EE 


jv=om=o + q2 + ( 2iV + i) 2 ^ 2 ]!-^ + q 2 + (2 M + 1) 2 7r 2 ] ‘ 
The non-zero coefficients of the double trigonometric series (IB44I) are: 

r 1/2 cl /2 

A NM {q) = dz dz' V(q,z,z r ) cos(2Nnz) cos(2AInz ), 

J- 1/2 J- 1/2 

q 2 F Snm 


(B44) 


8 q 5 F(l — coshg) cos(Nn) cos(M7r) 

(q 2 + 4N 2 n 2 ) 2 ' (q + smhq)(q 2 + 4N 2 ir 2 ) 2 (q 2 + 4M 2 n 2 ) 2 


(B45) 


+ 


and 


r L r* r 1 / 2 

BNAi(q) = / dz dz 1 V(q 1 z,z') sin[(2./V + l)7rz] sin[(2M + l)7r2:], 
J- 1/2 J- 1/2 

8 q 5 F(l + coshq) cos(Nn) cos(Mn) 


q 2 F S 


- 1/2 

NM 


(B46) 


(q 2 + (2N + 1) 2 7 t 2 ) 2 (g — sinhg)[g 2 + ( 2N + l) 2 7r 2 ] 2 [g 2 + (2 M + l) 2 7r 2 ] 2 


































